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The matter power spectrum at comoving scales of (1 — 40) /i~^Mpc is very sensitive to the presence 
of Warm Dark Matter (WDM) particles with large free streaming lengths. We present constraints 
on the mass of WDM particles from a combined analysis of the matter power spectrum inferred from 
the large samples of high resolution high signal-to-noise Lyman-a forest data of Kim et al. (2004) 
and Croft et al. (2002) and the cosmic microwave background data of WMAP. We obtain a lower 
limit of mwDM ~ 550 eV (2cr) for early decoupled thermal relics and mwDM ^ 2.0 keV (2cr) for 
sterile neutrinos. We also investigate the case where in addition to cold dark matter a light thermal 
gravitino with fixed effective temperature contributes significantly to the matter density. In that 
case the gravitino density is proportional to its mass, and we find an upper limit 7713/2 ~ 16 eV (2(t). 
This translates into a bound on the scale of supersymmetry breaking, Asusy < 260 TeV, for models 
of supersymmetric gauge mediation in which the gravitino is the lightest supersymmetric particle. 

PACS numbers: 98.80.Cq 
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I. INTRODUCTION 

It is now obscrvationally well established that the Uni- 
verse is close to flat and that matter accounts for about 
(25-30) percent of the total energy density most of which 
must be in the form of dark matter (DM) particles. Can- 
didates of dark matter particles are generally classified 
according to their velocity dispersion which defines a 
free-streaming length. On scales smaller than the free- 
streaming length, fluctuations in the dark matter density 
are erased and gravitational clustering is suppressed. The 
velocity dispersion of Cold Dark Matter (CDM) parti- 
cles is by definition so small that the corresponding free- 
streaming length is irrelevant for cosmological structure 
formation. That of Hot Dark Matter (HDM), e.g. Hght 
neutrinos, is only one or two orders of magnitude smaller 
than the speed of light, and smoothes out fluctuations 
in the total matter density even on galaxy cluster scales, 
which leads to strong bounds on their mass and density 
Q, 1^ Between these two limits, there exists an in- 
termediate range of dark matter candidates generically 
called Warm Dark Matter (WDM). If such particles are 
initially in thermal equilibrium, they decouple well be- 
fore ordinary neutrinos. As a result their temperature is 
smaller and their free-streaming length shorter than that 
of ordinary neutrinos. For instance, thermal relics with 
a mass of order 1 keV and a density VL^ ~ 0.25 would 
have a free-streaming length comparable to galaxy scales 
(Afs ~ 0.3 Mpc). 

There exist many WDM candidates whose origin is 
firmly rooted in particle physics. A prominent exam- 
ple is the gravitino, the supersymmetric partner of the 
graviton. The gravitino mass 7713/2 is generically of the 
order of Asusy /-^pi where Agusy is the scale of supersym- 
metry breaking and Mp is the Planckian scale. If Agusy is 



large the Lightest Supersymmetric Particle (LSP), which 
is then not the gravitino, decouples in the non-relativistic 
regime and provides a viable CDM candidate. If, how- 
ever, Asusy ^ 10^ GeV, as predicted by theories where su- 
persymmetry breaking is mediated by gauge interactions, 
the gravitino is likely to be the LSP. Such a light grav- 
itino has a wide range of possible masses (from lO^^eV 
up to the keV region). The gravitino - or better to say 
its helicity 1/2 component - decouples when it is still rel- 
ativistic. At this time the number of degrees of freedom 
in the Universe is typically of order one hundred, much 
larger than at neutrino decoupling. The effective grav- 
itino temperature is therefore always smaller than the 
neutrino temperature, and such light gravitinos can play 
the role of WDM Q . Their velocity dispersion is non- 
negligible during the time of structure formation, but 
smaller than the velocity dispersion of active neutrinos 
with the same mass. 

Other particles may decouple even earlier from thermal 
equilibrium while still rclativistic and may act as warm 
dark matter. The same is true for right-handed or ster- 
ile neutrinos added to the standard electroweak theory. 
Since their only direct coupling is to left-handed or ac- 
tive neutrinos, the most efficient production mechanism 
is via neutrino oscillations. If the production rate is al- 
ways less than the expansion rate, then these neutrinos 
will never be in thermal equilibrium. However, enough 
of them may be produced to account for the observed 
matter density. Right-handed (sterile) neutrinos with a 
mass of order keV are therefore natural WDM candidates 

The suppression of structures on scales of a Mpc and 
smaller, in WDM models with a particle mass of ~ 1 
kcV, has been invoked as a possible solution to two ma- 
jor apparent conflicts between CDM models and obser- 
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vations in the local Universe (see 0, H, S Q and ref- 
erences therein). The first problem concerns the inner 
mass density profiles of simulated dark matter halos that 
are generally more cuspy than inferred from the rota- 
tion curves of many dwarfs and low surface brightness 
galaxies. Secondly, N-body simulations of CDM models 
predict large numbers of low mass halos greatly in excess 
of the observed number of satellite galaxies in the Lo- 
cal Group. It appears, however, difficult to find WDM 
parameters which solve all apparent CDM problems si- 
multaneously and it is also not clear whether the tension 
between theory and observations is due to an inappro- 
priate comparison of numerical simulations and observa- 
tional data. The main difhculty is that the modeling of 
the inner density profile of galaxies and the abundance of 
satellite galaxies requires simulations on scales where the 
matter density field is highly non-linear at low rcdshift. 
Furthermore, many of the relevant baryonic processes are 
poorly understood. 

At z (2 — 3) the gravitational clustering of the mat- 
ter distribution spectrum at scales larger than the free 
streaming length of a thermal 1 keV WDM particle is still 
in the mildly non-linear regime. At these redshift high- 
resolution Lyman-a forest data offer an excellent probe 
of the matter power spectrum at these scales. Narayanan 
et al. [lol | compared the flux power spectrum and flux 
probability distribution of a small sample of high resolu- 
tion quasar (QSO) absorption spectra to that obtained 
from numerical dark matter simulations. In this way they 
were able to constrain the mass of a thermal WDM par- 
ticle (assumed to account for all the dark matter in the 
Universe) to be towdm ~ 750 cV. 

We will use here the linear DM power spectrum in- 
ferred from two large samples of QSO absorption spectra 
[iTl using state-of-the-art hydrodynamical simula- 
tions |l3l | combined with cosmic microwave background 
data from WMAP in order to obtain new constraints on 
the mass of possible thermally distributed WDM parti- 
cles. We thereby assume that WDM particles account for 
all the DM in our Universe and we do not put any priors 
on their temperature. We will also further investigate 
constraints on the mass of light gravitinos with realistic 
effective temperature. 

The paper is organized as follows. In Section II we will 
describe the effect of WDM on the power spectrum of the 
matter distribution. In Section III wc will summarize the 
method used in inferring the linear matter power spec- 
trum from Lyman-a forest observations and discuss the 
systematic errors involved. In Section IVA we will derive 
a lower limit for the mass of the dark matter particle in 
a pure A WDM model. In Section IVB, we will obtain an 
upper limit for the mass of light gravitinos in a mixed 
ACWDM universe. Section V contains our conclusions 
and some comments on the implications of our results 
for particle physics. 



II. WDM SIGNATURES ON THE MATTER 
POWER SPECTRUM 

Standard neutrinos arc believed to decouple from the 
plasma when they are still relativistic, slightly before 
electron-positron annihilation. They will keep the equi- 
librium distribution of a massless fermion with temper- 
ature ~ {A/nf/^T^. Their total mass summed over 
all neutrino species and their density are related through 
the well-known relation Q,^h^ ~ (TOi,/94eV) in the in- 
stantaneous decoupling limit. This result does, however, 
not apply to WDM for which: (i) either decoupling takes 
place earlier, (ii) or the particles were never in thermal 
equilibrium. 

If decoupling takes place at temperature To while 
there are g^iTjj) degrees of freedom the particles do not 
share the entropy release from the successive annihila- 
tions and their temperature is suppressed today by a 
factor 

/ 10.75 ^ 

Tu \9*{TDjJ 

as is the case e.g. for gravitinos. If the particles have 
never been in thermal equilibrium, there is a wide range 
of possible models. We will restrict our analysis here 
to the case of sterile neutrinos created from oscillations 
with active neutrinos, which share the same phase-space 
distribution but with an overall suppression factor x- In 
both cases, the number density of WDM is smaller than 
that of standard neutrinos, and the same present-day 
non-relativistic energy density corresponds to a larger 
particle mass, 

with (3 = {Tx/T^Y early decoupled thermal relics and 
P = X for sterile neutrinos. Due to their free-streaming 
velocity, WDM particles slow down the growth of struc- 
ture and suppress the total matter power spectrum P{k) 
on scales smaller than their free-streaming scale which is 
roughly given by 




see e.g. Note, however, that there exist different 

definitions of the free-streaming or damping scale. The 
effect of the free streaming on the matter distribution can 
be described by a relative "transfer function" , 

T{k) = [P(fc)AWDM/P(fc)ACDM]'/' (4) 

which is simply the square root of the ratio of the matter 
power spectrum in the presence of WDM to that in the 
presence of purely cold dark matter, for fixed cosmologi- 
cal parameters. 

Since we allowed for an arbitrary normalization factor 
X, the three parameters ojx, rnx and are independent. 
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Note, however, that the evolution equations for WDM as 
well as the corresponding terms in the Einstein equations 
can be reparametrized entirely in terms of two parame- 
ters, for instance uj^ and rrix/Tx (where Tx is expressed in 
units of the neutrino or photon temperature). There is a 
simple one-to-one correspondence between the masses of 
thermal WDM particles and sterile neutrinos for which 
the effect on the matter distribution and thus the transfer 
function for both models are identical |^. Writing equa- 
tion in the two cases and equating uj^ and m^jT^^ 
one obtains 



'Tistciiic 1. = 4.43 kcV 



/mthermaA 4/3/ 0.25(0.7)2 



V IkcV / 



(5) 

In the following, we will thus consider only the thermal 
model and simply translate any mass limit into that for 
the sterile neutrinos using eq. Q. 

Let us now discuss briefly some simple scenarios. 

Pure warm dark matter models. In the case in which 
the Universe contains only WDM (in addition to the 
usual baryon, radiation and cosmological constant com- 
ponents), the transfer function can be approximated by 
the fitting function 



r(fc) = [l + (afc)2'']-5/" 



(6) 



where a, the scale of the break, is a function of the 
WDM parameters, while the index v is fixed. With a 
Boltzmann code simulation (using either CMBFAST 
or CAMB 0), we find that ~ 1.12 is the best fit for 
k < 5 hMpc~^ and we obtain 



a = 0.24 



mx/J-x \ I i^x 



IkcV/T^y V 0-25(0.7)2 

/ irix \ 
VlkeV/ 



-0.16 

Mpc 



where the second line applies only to the case of thermal 
relics. This expression is close to that of (and exactly 
identical to the results of ^3 1 except for the front factor 
which was misprinted in this reference). 

In Figure ^ we plot the results of the numerical sim- 
ulations and compare them with our analytical fit. The 
power spectrum of the AWDM model differs from that 
of the corresponding ACDM model only at small scales. 
This is the reason why WDM has been suggested as a so- 
lution to the apparent "crisis" of CDM models on small 
scales. The matter density is fixed by the CMB and other 
data to LOx ~ 0.12. For pure AWDM models for which 
rux ~ 1 keV this gives g^iTu) ^10^. The particles there- 
fore decouple extremely early in such a model. Note that 
the light gravitino suggested by gauge mediated models 
of supersymmetry breaking decoupled when was much 
smaller, of the order of one hundred. 

Mixed models with CDM and gravitinos (or any ther- 
mal relics with g^(To) of a hundred). We will now con- 



sider the case of thermal relics which decoupled at a tem- 
perature when gt:(Tr)) ~ 100, as e.g. the light gravitino, 
which is likely to be the LSP in gauge-mediated theories 
of supersymmetry breaking. Note that g*{T]j) ^ 100 
corresponds to a decoupling temperature between a few 
GeV and a few TeV (possibly much higher, depending 
on which theory is assumed beyond the standard model). 
Note further that the temperature of thermal relics with 
g*{TD) ^ 100 is about a factor two smaller than that 
of ordinary neutrinos. They are thus not equivalent to 
HDM even in the small mass limit. For fixed g*, both 
the density ujx and the free-streaming wavenumber are 
proportional to the mass rUx. Taking g^ = 100, one gets 



0.117- 



rUx 



100 eV 



^FS 



1.5/iMpc"^ 



100 eV 



(8) 



In the absence of any cold dark matter, the preferred 
value ujx ~ 0.12 would lead to a mass of order 100 eV 
and to a complete smoothing out of fluctuations for 
k > 1.5/iMpc~^, at odds with the Lyman-a forest data 
(and with various other constrains related to structure 
formation 0|). However, in the presence of CDM, 
small-scale structures are not completely suppressed and 
the free-streaming effect leads to a step-like transfer func- 
tion (as for massive neutrinos, but with a step at a 
much smaller scale), while in the small limit stan- 
dard ACDM is recovered. For a light gravitino in the 
the mass range 10 eV ^ ^3/2 ~ 100 eV, 5* (To) varies 
between 87 and 101 depending on the details of the su- 
persymmetric model [13- As we already pointed out a 
light gravitino with larger mass would result in a matter 
density which exceeds that observed and as we will show 
later smaller masses are hard to constrain with current 
observations. In Figure 1 we show the transfer function 
computed numerically for a few values of the mass in the 
range (10 — 100) eV, corresponding to fcpg in the range 
(0.15 — 1.5)/iMpc~^. Since the matter power spectrum 
measured from galaxy redshift surveys can be compared 
to linear predictions up to fc ~ 0.15 /iMpc, it has little 
sensitivity to this model but as we will show in Section 
IV the matter power spectrum inferred from Lyman-a 
data provides a tight upper limit on the gravitino mass. 

Mixed models with warmer thermal relics. Finally, the 
WDM particles could decouple at temperature of order 
1 GeV or below (i.e. just before or after the QCD phase 
transition), when the number of degrees of freedom was 
in the range 10.75 < < 80 leading later to a tem- 
perature 0.5 < {Tx/Tu) < 1. This case smoothly inter- 
polates between that of ordinary neutrinos (HDM) and 
that of the previous paragraph. The free-streaming scale 
is smaller than for the gravitino and these models have 
been constrained by the matter power spectrum recon- 
structed from galaxy redshift surveys 0]- We will not 
consider this case here. 
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FIG. 1: The WDM transfer function T{k) defined in 
equation (|1J for various models (all with the same cos- 
mological parameters Qb = 0.05, Odm ~ 0.25, f^A = 
0.70). (a) Pure AWDM model with, from left to right, 
Trc/T^ = 0.5,0.3,0.25,0.226,0.2, corresponding to = 
92, 427, 738, 1000, 1441 eV. The solid curves are obtained nu- 
merically, the long-dashed curves (essentially indistinguish- 
able from the solid curves) show the analytical fits based 
on equation O, the short-dashed curves those based on 
Ref. PI- (b) Mixed AC WDM model for a warm compo- 
nent which decoupled when gt{TD) = 100, like e.g. a light 
gravitino (with mass proportional to density and f2cDM = 
0.25 — i}x). The solid curves show the numerical results for 
= 10, 20, 30, 50, 70, 100 eV, from top right to bottom left. 
At the top-axis we show the wavenumber scale in s/km (as- 
suming the above cosmological model and z = 2.5 interme- 
diate between the two Lyman-o? data sets analysed here). In 
both panels the double arrow indicates the range of wavenum- 
bers used in our analysis. 



III. PROBING THE MATTER POWER 
SPECTRUM WITH THE LYMAN-a FOREST IN 
QSO ABSORPTION SPECTRA 

It is well established by analytical calculation and hy- 
drodynamical simulations that the Lyman-a forest blue- 
ward of the Lyman-a emission line in QSO spectra is 
produced by the inhomogenous distribution of a warm 
(~ 10'' K) and photoionized intergalactic medium (IGM) 
along the line of sight. The opacity fluctuations in the 
spectra arise from fluctuations in the matter density and 
trace the gravitational clustering of the matter distribu- 
tion in the quasi- linear regime [I3 ■ The Lyman-a forest 
has thus been used extensively as a probe of the matter 
power s pec trum on comoving scales of (1 — 40) h~^Mpc 

The Lyman-a optical depth in velocity space u (km/s) 
is related to the neutral hydrogen distribution in real 



t{u) 



Hiz) 



dy nmiy) V [u - y - vii {y) , b{y)] dy 



(9) 

where a^ ^ — 4.45 x 10 cm^ is the hydrogen Lya cross- 
section, y is the real-space coordinate (in km s~^), V 
is the standard Voigt profile normalized in real-space, 
b = (2fcBT/mc2)i/2 is the velocity dispersion in units of 
c, H{z) the Hubble parameter, nni is the local density 
of neutral hydrogen and v\\ is the peculiar velocity along 
the line-of-sight. The density of neutral hydrogen can 
be obtained by solving the photoionization equilibrium 
equation (e.g. 0|). The neutral hydrogen in the IGM 
responsible for the Lyman-a forest absorptions is highly 
ionized due to the metagalactic ultraviolet (UV) back- 
ground radiation produced by stars and QSOs at high 
redshift. This optically thin gas in photoionization equi- 
librium produces a Lyman-a optical depth of order unity. 

The balance between the photoionization heating by 
the UV background and adiabatic cooling by the expan- 
sion of the universe drives most of the gas with 5b < 10, 
which dominates the Lyman-a opacity, onto a power-law 
density relation T = Tq (1 + Sb)'^~^ , where the param- 
eters To and 7 depend on the reionization history and 
spectral shape of the UV background and Si, is the local 
gas overdensity {1 + Sb = Pb / Pb) ■ 

The relevant physical processes can be readily mod- 
elled in hydro-dynamical simulations. The physics of 
a photoionized IGM that traces the dark matter distri- 
bution is, however, sufficiently simple that considerable 
insight can be gained from analytical modeling of the 
IGM opacity based on the so called Fluctuating Gunn 
Peterson Approximation neglecting the effect of peculiar 
velocities and the thermal broadening The Fluc- 

tuating Gunn Peterson Approximation makes use of the 
power-law temperature density relation and describes the 
relation between Lyman-a opacity and gas density (see 
10,13) along a given line of sight as follows. 



t(z) oc (l + Sb{z)fT-''-' 
A{z) = 0.433 ^ 



(2) = Aiz){l + 6b{z)f , (10) 



0.65 



3.5 

-1 



nhh^ 



0.02 

H{z)/Ho 
3.68 



f To 



V6000K 



Fhi 



1.5 X 10-12 s-i 



where /3 = 2-0.7 (7-I) in the range 1.6-1.8, Fhi the HI 
photoionization rate, Hq = h 100 km/s/Mpc the Hubble 
parameter at redshift zero. For a quantitative anlysis, 
however, full hydro-dynamical simulations, which prop- 
erly simulate the non-linear evolution of the IGM and its 
thermal state, are needed. 

Equations H9I10|I show how the observed flux F = 
exp (— r) depends on the underlying local gas density pb, 
which in turn is simply related to the dark matter den- 
sity, at least at large scales where the baryonic pressure 
can be neglected |23l- Statistical properties of the flux 
distribution, such as the flux power spectrum, are thus 



5 



closely related to the statistical properties of the under- 
lying matter density field. 



A. The data: from the quasar spectra to the flux 
power spectrum 

The power spectrum of the observed flux in high- 
resolution Lyman-a forest data provides meaningful con- 
straints on the dark matter power spectrum on scales 
of 0.003 s/km < k < 0.03 s/km, roughly corresponding 
to scales of (1 — 40) /i^^Mpc (somewhat dependent on 
the cosmological model). At larger scales the errors due 
to uncertainties in fitting a continuum (i.e. in removing 
the long wavelength dependence of the spectrum emit- 
ted by each QSO) become very large while at smaller 
scales the contribution of metal absorption systems be- 
comes dominant (see e.g. ^3,|2^). In this paper, we will 
use the dark matter power spectrum that Viel, Haehnelt 
& Springel [H (VHS) inferred fro m the flux power spec- 
tra of the Croft et al. [13 (C02) sample and the LUQAS 
sample of high-resolution Lyman-a forest data . The 
C02 sample consists of 30 Keck high resolution HIRES 
spectra and 23 Keck low resolution LRIS spectra and 
has a median redshift of z = 2.72. The LUQAS sample 
contains 27 spectra taken with the UVES spectrograph 
and has a median redshift of 2; = 2.125. The resolution 
of the spectra is 6 km/s, 8 km/s and 130 km/s for the 
UVES, HIRES and LRIS spectra, respectively. The S/N 
per resolution element is typically 30-50. Damped and 
sub-damped Lyman-a systems have been removed from 
the LUQAS sample and their im pac t on the flux power 
spectrum has been quantified by |l2l |. Estimates for the 
errors introduced by continuum fitting, the presence of 
metal lines in the forest region and strong absorptions 
systems have also been made 0, 0, |2^ • 

We stress that the total sample that went into the es- 
timate of the linear power spectrum used here is about 
10 times larger than the one used in the study by p^. 



B. From the flux power spectrum to the linear 
matter power spectrum 

VHS have used numerical simulation to calibrate the 
relation between fiux power spectrum and linear dark 
matter power spectrum with a method proposed by C02 
and improved by and VHS. A set of hydro-dynamical 
simulations for a coarse grid of the relevant parameters is 
used to find a model that provides a reasonable but not 
exact fit to the observed fiux power spectrum. It is then 
assumed that the differences between the model and the 
observed linear power spectrum depend linearly on the 
matter power spectrum. 

The hydro-dynamical simulations are used to define a 
bias function between fiux and matter power spectrum: 
Ppik) = 6^(fc) P{k), on the range of scales of interest. 
In this way the linear matter power spectrum can be re- 



covered with reasonable computational resources. This 
method has been found to be robust provided the system- 
atic uncertainties are properly taken into account [l^l29l |. 
Running hydro-dynamical simulations for a fine grid of 
all the relevant parameters is unfortunately computation- 
ally prohibitive. 

The use of state-of-the-art hydro-dynamical simula- 
tions is a significant improvement compared to previous 
studies which used numerical simulation of dark mat- 
ter only [lol [1^ . We use a new version of the paral- 
lel TreeSPH code GADGET ^ in its TreePM mode 
which speeds up the calculation of long-range gravita- 
tional forces considerably. The simulations are performed 
with periodic boundary conditions with an equal number 
of dark matter and gas particles. Radiative cooling and 
heating processes are followed using an implementation 
similar to [23 for a primordial mix of hydrogen and he- 
lium. The UV background is given by [31|. In order to 
maximise the speed of the simulation a simplified crite- 
rion of star formation has been applied: all the gas at 
overdensities larger than 1000 times the mean overden- 
sity is turned into stars [l^ . The simulations were run on 
COSMOS, a 152 Gb shared memory Altix 3700 with 152 
CPUs hosted at the Department of Applied Mathematics 
and Theoretical Physics (Cambridge). 

In order to check the impact of a WDM particle on the 
flux power spectrum, we have run a set of simulations 
with different WDM particle masses and compared them 
with a ACDM simulation with the same phases of the 
initial conditions. The cosmological parameters for both 
simulations are Om = 0.26, SIa = 0.74, fie = 0.0463 and 
Ho = 72kms-iMpc-^ The ACDM transfer functions 
have been computed with CMBFAST [l^. The transfer 
function for the WDM model has been adopted from |3l . 
The simulations were run with a box size of 30 comoving 
h^^ Mpc with 2 X 200^^ gas and dark matter particles. 

In Figure m we show the fractional difference between 
the flux power spectra of a AWDM model and that of 
a ACDM model for a range of WDM particle masses at 
z = 2.5, 3, 4. The free-streaming of the WDM particles 
has little effect for particles with masses ^ 1 keV (ther- 
mal WDM) at scales for which the linear matter power 
spectrum can be reliably inferred (less than 5% at z=2.5 
increasing to 15% at z=4). For smaller masses of the 
WDM particle there is. however, a dramatic reduction of 
the ffux power spectrum. Note that at small small scales 
^ 0.1 s/km, the ffux power spectrum shows a bump. This 
is the effect of the non-linear evolution of the matter dis- 
tribution which results in a transfer of power from large 
scales to small scales |l3 . The diamonds in Figure|21show 
the effect of adding randomly oriented velocities drawn 
from a Fermi-Dirac distribution in the initial conditions 
for the model with mwDM=l keV. We confirm the result 
of ^3 that the effect of adding these free-streaming ve- 
locities on the fiux power spectrum is small, of the order 
of 2 % for k < 0.03 s/km. We will neglect the effect in 
the following. 

Figures n ^nd 121 suggest that the Lyman-a forest data 
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FIG. 2: Fractional difFerence of the flux power spectrum for a hydro-dynamical simulation of a AWDM model and the flux 
power spectrum for a ACDM model. The simulations have a box size of 30 Mpc and 200'^ gas and 200'^ dark matter 
particles. The results are for three different redshifts [z — 4, 3, 2.5, respectively from left to right) and for five different values 
of the (thermal WDM) particle mass: 0.5 keV (thick dotted), 1 keV (thin dashed), 2 keV (thin dotted), 4 keV (thick dashed) 
and 8 keV (continuous). The diamonds show the results for a simulation with a WDM particle of mass IkeV when random 
oriented velocity dispersion drawn from a Fermi-Dirac distribution have been added in the initial conditions. 



will provide a tight lower limit to the mass of WDM 
particles. We will quantify this limit in the next sec- 
tion using the linear dark matter power spectrum in the 
range (0.003—0.03) s/km as inferred by VHS. We thereby 
checked how strongly the bias function used to infer the 
linear matter power spectrum depends on the presence of 
WDM. For this purpose we have run a simulation with 
a. 60 comoving Mpc box with 2 x 400^ gas and wark 
dark matter particles for a (thermal WDM) particle mass 
"T-WDM = 1 keV. In the relevant wavenumber range, the 
difference between the bias function extracted from this 
simulation and that of a ACDM model (B2 model in 
VHS) is less than 3%, as we can see from Figure O For 
smaller WDM particle masses the difference is expected 
to become larger but as we discuss in Section FlV Al this 
will have little impact on our final results. 



C. Systematics Errors 

There is a number of systematic uncertainties and sta- 
tistical errors which affect the inferred power spectrum 
and an extensive discussion can be found in [T3 . [l3L[2^ . 
VHS estimated the uncertainty of the overall rms fluctu- 
ation amplitude of matter fluctuation to be 14.5 % with 
a wide range of different factors contributing. 

In the following we present a brief summary. The ef- 
fective optical depth, TeS = — In < F > which is essen- 
tial for the calibration procedure has to be determined 
separately from the absorption spectra. As discussed in 




0.05 . I 1 

0.07 0.70 

k (s/km) 

FIG. 3: Bias function 6^(A:) between the flux power spec- 
trum and the linear dark matter power spectrum, Ppik) = 
b^{k) P{k), for two cosmological models: ACDM (diamonds) 
and AWDM (squares), with a (thermal WDM) particle mass 
jnwDM = 1 keV. The simulation has a boxsize of 60 comov- 
ing h~^Mpc and 2 x 400'' (gas and dark matter) particles (see 
Section fill Bl for further details on other simulation parame- 
ters). The shaded area indicates the range of wavenumbers 
used in our analysis. 
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VHS; there is a considerable spread in the measurement 
of the effective optical depth in the literature. Deter- 
minations from low-resolution low S/N spectra give sys- 
tematically higher values than high-resolution high S/N 
spectra. However, there is little doubt that the lower val- 
ues from high-resolution high S/N spectra are appropri- 
ate and the range suggested in VHS leads to a 8% uncer- 
tainty in the rms fluctuation amplitude of the matter den- 
sity field (Table 5 in VHS). Other uncertainties are the 
slope and normalization of the temperature-density rela- 
tion of the absorbing gas which is usually parametrised 
as T = To (1 -|- 6b)'^~'^. To and 7 together contribute up 
to 5% to the error of the inferred fluctuation amplitude. 
VHS further estimated that uncertainties due to the C02 
method (due to fitting the observed flux power spectrum 
with a bias function which is extracted at a slightly differ- 
ent redshift than the observations) contribute about 5%. 
They further assigned a 5 % uncertainty to the somewhat 
uncertain effect of galactic winds and finally an 8% uncer- 
tainty due the numerical simulations (codes used by dif- 
ferent groups give somewhat different results). Summed 
in quadrature, all these errors led to the estimate of the 
overall uncertainty of 14.5% in the rms fluctuation am- 
plitude of the matter density field. 

For our analysis we use the inferred DM power spec- 
trum in the range 0.003 s/km < k < 0.03 s/km as given 
in Table 4 of VHS. Note that we have reduced the values 
by 7% to mimick a model with 7 = 1.3, the middle of the 
plausible range for 7 [s^l . Unfortunately at smaller scales 
the systematic errors become prohibitively large mainly 
due to the large contribution of metal absorption lines to 
the flux power spectrum [T^ (their Figure 3) and due to 
the much larger sensitivity of the flux power spectrum to 
the thermal state of the gas at these scales. 



IV. BOUNDS ON THE WDM MASS FROM 
WMAP AND THE LYMAN-q FOREST 

In this section we will show the results of our com- 
bined analysis of CMB data and Lyman-a forest data. 
The CMB power spectrum has been measured by the 
WMAP team over a large range of multipoles {I < 800) to 
unprecedented precision on the full sky [l| . For our anal- 
ysis we use the first year data release of the WMAP tem- 
perature and temperature-polarization cross-correlation 
power spectrum. 

For the range of parameters nix and lu^ considered 
here, CMB anisotropics are not sensitive to the veloci- 
ties of WDM particles and can therefore not discriminate 
between WDM and CDM models. It is nevertheless cru- 
cial to combine the CMB and Lyman-a data in order to 
get some handle on other cosmological parameters. Since 
all recent CMB experiments point toward a flat Universe 
with adiabatic scalar primordial fluctuations, we adopt 
these results as a prior. We treat the three ordinary 
neutrinos as massless (including small masses would not 
change significantly our results). We use the publicly 



available code CAMB [l^ in order to compute the theo- 
retical prediction for the C; coefficients of the CMB tem- 
perature and polarization power spectra, as well as the 
matter spectrum P{k). We derive the marginalized like- 
lihood of each cosmological parameter with Monte Carlo 
Markov Chains as implemented in the publicly available 
code CosmoMC js^- The calculations were performed 
on the Cambridge COSMOS supercomputer. We extended 
the CosmoMC package in order to include the Lyman- 
a forest data in a way similar to Ref. [s^l ■ 

A. Lower bounds on rUx in a pure AWDM 
Universe. 

We first compute numerically the ACDM matter power 
spectrum, and then use equation (O with = 1.2 to ob- 
tain the pure AWDM matter power spectrum leaving the 
break scale a as a free parameter. Apart from a, our pa- 
rameter set consists of the six usual parameters of the 
minimal cosmological model: the normalization and tilt 
of the primordial spectrum, the baryon density, the dark 
matter density wdm (equal here to uJx), the ratio of the 
sound horizon to the angular diameter distance, and the 
optical depth to reionization. Note that the Hubble pa- 
rameter and the cosmological constant can be derived 
from these quantities, as well as as (the parameter de- 
scribing the amplitude of matter fluctuations around the 
scale R = 8/i~^Mpc). In addition to these seven cosmo- 
logical parameters (with flat priors), we vary the param- 
eter A (see 0,0) which accounts for the uncertainty 
of the overall amplitude of the matter power spectrum 
inferred from the Lyman-a forest data (with a Gaussian 
prior as described in H^l). 





AWDM 


ACWDM 




0.124 ±0.015 


0.149 ±0.019 




0.024 ± 0.001 


0.024 ± 0.001 


h 


0.72 ± 0.06 


0.71 ±0.06 


T 


0.18 ±0.09 


0.17 ±0.08 


0-8 


0.96 ±0.08 


0.86 ± 0.09 


n 


1.01 ±0.04 


1.00 ±0.04 


a (Mpc/h) 


0.06 ± 0.03 








0.05 ± 0.04 



TABLE I: The marginalized results for the recovered cosmo- 
logical parameters for our WMAP + Lyman-a runs. The left 
column is for pure AWDM models while the second column is 
for mixed ACWDM models (with a fixed WDM temperature 
corresponding to (Td) = 100). The values correspond to 
the peaks of the marginalized probability distributions. Er- 
rors are the 68% confidence limits. 

The best-fit values for the cosmological parameters are 
summarized in the first column of Table 1. In Figure 01 
(first panel), we show the marginalized likelihood for the 
break scale a. The probability peaks at 0.07 /i~^Mpc, 
but the preference for a non-vanishing a is not statisti- 
cally significant, and our main result is the upper bound 
a ^ 0.11/i~^Mpc at the 2a confidence level. Using 
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eq. (O, we can derive the likelihood for the mass rrix, or 
better, for its inverse, in order to keep a nearly flat prior. 
The second panel of Figure 0] shows the two-dimensional 
likelihoods for the WDM parameters {m~^ , uj^), assum- 
ing the WDM is a thermal relic. In this case the 2a 
bound on the mass is mx ^ 550 eV, while for a non- 
thermal sterile neutrino as described in Section II this 
limit translates into rrix <i 2.0 kcV. 

Note that for particle masses ^ 0.5 keV the difference 
of the linear matter power spectrum of a WDM model 
compared to that of the corresponding CDM model are 
not small anymore at the relevant scales (Fig. 1). This 
makes the use of the linear bias relation as described in 
Section III questionable when deriving the likelihoods. 
The sensitivity to the mass of the WDM particle is how- 
ever so strong that this should have no effect on our con- 
clusions other than a possible small shift of our 1 and 2 
a bounds in this region. 

We find no strong correlation between the mass of the 
WDM particle and other cosmological parameters. This 
is expected as the characteristic break due to the free- 
streaming length of the WDM particles cannot be easily 
mimicked by the effect of other parameters. Note how- 
ever that the constraints on flxh'^, cr^ and the tilt Us of 
the primordial spectrum are significantly weakened for 
WDM particle masses ^ IkeV and that somewhat larger 
values of cts and Us become allowed. 



B. Upper bound on irix in a ACWDM Universe 




FIG. 4: Results for the pure AWDM model. The top left 
panel shows the likelihood for the break scale a. The other 
three panels show the la and 2a contours for (m^^, lox), 
(m~^, as) and (m~^, Ub), repectively, where nix is the mass 
of a thermal WDM particle. The solid curves correspond to 
the full marginalized likelihoods, while the dashed curves and 
color/shade levels show the mean likelihood of the samples in 
the Markov chains. 



Let us now perform a similar analysis for the second 
model described in Section II: a ACWDM Universe with 
a warm thermal relic assumed to decouple when the num- 
ber of degrees of freedom was g^, ~ 100 (as for instance 
for a light gravitino playing the role of the LSP). We 
choose the same set of cosmological parameters as in 
the previous subsection, except that the WDM sector is 
not described by the scale a but by the density fraction 
jx = W2;/wdm (where wdm = + wcdm)- The mass rrix 
is proportional to the density lOx (Eq- lU ; a-nd the limit 
rUx — > corresponds to the standard ACDM model. In 
this case, we will therefore obtain an upper instead of a 
lower limit on the mass. 

In Figure ^ (first panel), we show the marginalized 
likelihood for the WDM fraction. The 2a Bayesian con- 
fidence limit ]x 0.12 can be translated into a limit 
on the mass, nix 16 eV. The other three panels show 
two-dimensional contours for and the most correlated 
variables: cts, lox and Ug. There is now some correlation 
between the mass and the total dark matter density, be- 
cause increasing t^DM tilts the small-scale matter power 
spectrum in such a way that it compensates to some ex- 
tent the red tilting introduced by mx- For the same rea- 
son, one could expect some correlation with the tilt of 
the primordial spectrum. This does not occur because 
the value of the tilt is fixed more precisely than that of 
wdm by the CMB data. The gravitino model predicts 



globally less power on small scales than ACDM. It there- 
fore prefers smaller values of (Tg in a combined analysis 
with the CMB data on large scales. Precise indepen- 
dent measurements of erg could thus at least in principle 
tighten the limit on the mass of a putative gravitino. 

Note that in this analysis, we assumed (7*(Td) = 100, 
while as mentioned in Section II detailed studies of real- 
istic supersymmetric models predict 87 < g*{Ti}) < 104, 
for the mass range investigated here 0|. We have re- 
peated our analysis assuming {To ) = 90 (which results 
in a slightly smaller coefficient of proportionality between 
the WDM mass and density) and found no significant dif- 
ference. 



V. CONCLUSIONS 

We have obtained constraints on the mass of ther- 
mal and non-thermal WDM particles using a combined 
analysis of the cosmic microwave background data from 
WMAP and the matter power spectrum inferred from a 
large sample of Lyman-a forest data. 

Assuming that the dark matter is entirely in the form 
of warm thermal relics, we find a 2a upper bound irix ^ 
550 eV on the mass of the DM particle. This confirms the 
result of Ref. ^3 who got a slightly bigger lower limit. 
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FIG. 5: Results for the ACWDM model with a light gravitino 
(or a particle decoupling when <;« = 100). The top left panel 
shows the likelihood for the WDM density fraction . Other 
panels show the two-dimensional contours for the gravitino 
mass ra^ and the most correlated variables: erg, and lox- 
The solid curves correspond to the full marginalized likeli- 
hoods, while the dashed curves and color/shade levels show 
the mean likelihood of the samples in the Markov chains. 



Wc have used a significantly larger sample of QSO spec- 
tra, applied an improved analysis which used state-of- 
the-art hydro-dynamical simulations and taken into ac- 
count a wide range of systematic uncertainties giving us 
confidence that our result is robust. Pushing this lower 
limit to larger masses will require accurate modeling of 
the Lyman-a forest data at smaller scales. This is diffi- 
cult due to the sensitivity of the flux power spectrum to 
metal absorption and the thermal history of the gas at 
these scales. 

For non-thermal WDM particles like e.g. a sterile neu- 
trino (with the same phase-space distribution as that of 
a standard neutrino, modulo a global suppression factor) 
this limit translates into m ^ 2.0 keV (2(7). Note that in 
this case there exists also an upper bound m ^ h keV, 
from limits on the sterile neutrino decay rate in dark 
matter halos set by X-ray observations [33 ■ It has been 
noticed that reionization at z ~ 6 requires mwDM^?^ 1 
keV, in addition if Swift discovers z ~ 10 gamma ray 
bursts then this would require even more stringent lower 
limits 36] . The existence of a WDM particle with a mass 
not much larger than these limits would increase the in- 
ferred amplitude of the matter power spectrum erg and 
the inferred tilt of the primordial spectrum compared to 
a CDM cosmology. 



Our results can be easily converted into a bound on 
a possible tiny interaction between CDM and photons, 
using the analysis of Ref. |33| (who showed that such an 
interaction would affect structure formation in the same 
way as WDM). 

We have also obtained constraints on the mass of ther- 
mal relics which decoupled at temperatures of the order 
of GeV or TeV and on their contribution to the total mat- 
ter density. In this case the number of degrees of free- 
dom in the Universe at the time of decoupling (7*(Td) 
was of order of one hundred and the particles cannot 
make up all the dark matter in the Universe ~ otherwise 
they would have a mass of order 100 eV which would 
contradict the limit discussed above. A light gravitino 
is a prime example of such a particle and would only 
be viable if it coexisted with ordinary cold dark matter 
(or with another form of warm dark matter particle with 
significantly larger mass). In this case, for particles with 
9*{Td) in the range from 90 to 100, we find a 2a bound 
nix ^ 16 eV. So, if the gravitino is the LSP, as sug- 
gested by gauge-mediated susy breaking scenarios, the 
susy breaking scale is limited from above. 



A,, 



1/2 



< 260 TeV 



(11) 



This conclusion is robust and relevant for supersymmetry 
searching at future accelerator machines as LHC. Indeed, 
if we assume that after decoupling the gravitino back- 
ground was enhanced by the decay of the NSP (Next -to- 
lightest Siipersymmetric Particle), as it happens in some 
scenarios |38j . the total density of gravitinos as well as 
their velocity dispersion can only be enhanced, making 
the bound even stronger. The only possible way to evade 
this bound on the susy breaking scale (still assuming that 
the gravitino is the LSP) would be to assume some en- 
tropy production after gravitino decoupling, as suggested 
in Refs. [MliS. 
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